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1 Introduction 



During the past several years, there has been significant progress in the determination of 
the values of the strong-coupling constant as at the Z-boson mass scale and the heavy- 
quark masses mh in perturbative quantum chromo-dynamics (QCD). The matching rela- 
tions which determine the finite discontinuities of at the heavy-quark thresholds in the 
modified minimal-subtraction (MS) scheme and are necessary for a precise determination 
of Us from a global data fit were recently obtained in semi- analytical [1, and analytical 
[2| form at the four-loop level. Furthermore, there is significant progress in multi-loop 
technology aiming towards an analytical evaluation of the QCD j3 function at five loops. 

In order to precisely determine heavy-quark masses with the help of QCD sum rules 
0, it is necessary to know in detail the heavy-quark contribution to the QCD vacuum 
polarisation function. The most important ingredients for this analysis are the lowest 
moments of the polarisation function (see Ref. Indeed, it was found that the first 
moment is best suited for this analysis, since it has the weakest dependence on non- 
perturbative effects and the details of the threshold region, thus leading to the smallest 
theoretical uncertainty [5]. 

At the four-loop level, the study of the first two moments has been started in Ref. ^ 
with calculations of diagrams including two internal loops of massive and massless quarks 
coupled to gluons. The results in Ref. contain two numerical constants, denoted as 
A^io and A'"2o, which were not expressed in terms of basic transcendental numbers. In the 
meantime, the corresponding Feynman diagrams have been calculated analytically |7|8j . 
The terms proportional to ainj~'^, where rii in the number of light quarks, are known in 
all orders j of perturbative QCD References |l()pilj contain the complete four-loop 
contributions to the first two Taylor coefficients from non-singlet diagrams. The singlet 
contributions were studied in Refs. |lH12j . 

The results of Refs. |imilj contain some four-loop tadpoles which are only given in 
numerical form. The purpose of this Letter is to evaluate them in terms of standard 
transcendental numbers so as to represent the first two moments of the heavy-quark 
vacuum polarisation function at four loops in completely analytical form. 

The content of this Letter is as follows. Section |21 contains the basic formulae. In 
Section El we present our analytic results for the tadpoles that, in Refs. |lU|13j . are 
called T54, Tq2, and T91 as well as for the first two moments of the heavy-quark vacuum 
polarisation at four loops in QCD. A brief summary is given in Section 0] 

2 Approach 

Since the vector current j^{x) = \E'(x)7''\E'(x) constructed from the heavy-quark field ^(x) 
is conserved, the correlator. 



U^''{q)=i j dxe"'^0\Tf{x)f{0)\0) 

= (-gV^ + gV)n(g'), 



(1) 
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can be expressed in terms of a single scalar function 11 (g^). The latter is of great phe- 
nomenological interest because it is related to the experimental observable 

„ , , cr(e+e^ hadrons) 

Ris) = \ ^ _ 

= 127rlmn(s + 2e), (2) 

where s is the square of the e^e" centre-of-mass energy. Equation Q is equivalent to an 
infinite number of equalities between experimental moments, 

= I ^^(^)' (3) 



and their theoretical counterparts 

where Ch is the electric-charge quantum number of the heavy quark, Wih = rrih {rrih), and 
Cn are the coefficients of the Taylor expansion 

^('l') = ^.j:CnZn (5) 

n>0 

mz = g^/(4m^). The perturbative calculation of 

where a,^ = ctg i^h), only involves tadpole diagrams and is, therefore, considerably simpler 
than the one of n(g^), which has a complicated dependence on 1. Here, it is understood 
that asi/J') and mh{^) are defined in the MS renormalisation scheme, which is based on 
dimensional regularisation with D = 4 — 2e space-time dimensions and 't Hooft mass scale 
/i. 

Thanks to the strong hierarchy among the quark masses, one may split the number 
Uf = rii + Uh of quark fiavours into rih = 1 massive ones and n/ massless ones. In the 
following, we keep the variable Uh generic. 



3 Results 

In our previous paper p], we introduced a technique that allows one to analytically 
evaluate a large class of four-loop tadpole diagrams with one non-zero mass. Specifically, 
the considered diagrams are transformed into integral representations whose integrands 
contain only one-loop tadpoles with new propagators having masses that depend on the 
variables of integration. This technique is based on the differential equation method jH] . 
A similar technique was applied to certain types of two-loop diagrams in Ref. [TK] . 
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Figure 1: Massive four- loop tadpole diagrams T54, Tq2, and T91. 



To illustrate the usefulness of this technique, in Ref. 0, we also evaluated the tadpoles 
that are denoted in Ref. [H] as A^'io and N20. The finite parts of our results agree with 
the numerical results of Ref. jOj and the analytical ones of Ref. [7], while our analytical 
results for the 0(e) terms are new. 

These results together with those from Ref. j2j offer us the opportunity to obtain 
analytical results for terms of the expansion in e of the master integrals denoted as T54, 
Tq2, and T91 in recent papers |lU|13j for which only numerical values have been available 
so far. The corresponding Feynman diagrams are depicted in Fig. Specifically, Tq2 can 
be extracted from A^'io jHJ using Eqs. (16) and (18) of Ref. jB]. This result can then be 
used in combination with T91 [2] to extract T54 from Eq. (50) of Ref. fT^. The results 
read 

1 9 415 4991 89383 / 3367679 1792^,^, 

+ M =;;;;7^ + ^C(3) 



N 2e3 36e2 216e 1296 V 7776 9 

of 137735095 43527r^ 68992^,^, 1024, 
+ + C(3) + h 

V 46656 135 27 ^ 3 

4908181487 1 675527r^ 174087r\ 1602496^,, 
+ + In 2 + C 3 

V 279936 405 45 81 ^ 

_5!!!5c(5) + 5^.. + 4096O+O(.'), (7) 
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.4 



+ e (-1156 - ^ + ^C(3) + 96C(5) + 6464 

2/ 20 9 3 8 71 07r^ 24167r\ Sn^ 12952^,, 

+ + n 2 + C 3 

V 3 9 45 21 3^^^ 

64 \ 
+ yC'(3) - 2400C(5) + 704^4 + 51265 j + 0(e=^), 
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iro_4 Q-TQ 

T91 = In 2 + — C(5) - 4865 + 0(e), (9) 

where 

h = -^(tt' - In^ 2) In^ 2 + 8 Lu Q ) , 

h = ^(57r2 - 3 In^ 2) ln=^ 2 + 8 Lig , (10) 

and = r'^(l + e)(m^)~^'^ is a normalisation factor. The choice = 1 corresponds to 
the normahsation adopted in Ref. JH]- The terms through O(e^) in Eq. ((7j) and those 
through 0(e) in Eq. (jH} agree with Eqs. (38) and (42) in Ref. [1^], respectively, while the 
other terms represent new results. 

Equations extend the set of four-loop tadpole master integrals with one non- 

vanishing mass known in fully analytical form, and represent essential ingredients for 
high-precision predictions of numerous observables of current phenomenological interest. 
After changing the overall normalisation, we find agreement with the numerical values 
given in Ref. 0. We note in passing that also the other master integrals presented in 
Ref. can be written in a more compact form if the combinations 64 and 65 of Eq. (jlUj) 
are introduced. 

Equation Q is the key result of Ref. j2j and a crucial ingredient for the four-loop 
matching condition of as{fi) at the heavy-quark thresholds |1|2I . Equations (0)-© are 

the missing links for a fully analytical representation of the moments C^q^ and C^-^^ |lUj . 
which is the main goal of this paper. 
Our results for and cf^ read 



-(3) 2 17897 31 \ ( 7043 497r^ 

127 1 , A . 2 f 610843 661 



34992 12960 

C 3 64 \+ni[ C 3 

324^^ ' 36 V ^ V 2449440 2835^^ ' 

( 7162 9 8 5 3371^ 21343,,^, 25, 
+ ni \ C(3) &4 

V 46656 116640 3888 ^^^ 324 

/ 83433703 14873?!^ 14509529^, , 5,, , 
+ Uh + C(3) + -C(5) 

V 8164800 18225 340200 ^ ^ 3^ ^ 

7091, \ 5266559 7029597r^ 12897r^ 

64 1 - — — — — ln2 



'5, 



810 V 466560 699840 3645 
1688407,,^, 89,,^, 313, 328, 
+ ^296r^(')-24^^')- 1944^^-^^ 



-(3) 2/^42173 112^, A ^ /262877 142l7r4 
^\ =^l\ 7^^7777 - 77T7C(3) + niiiu ^^^^^^ + 



V98415 405^' 7 V787320 174960 

38909 29 \ ^2 A63868 3287 ' 

58320^^ ~ 486 7 ^ V 295245 " 7290^^ \ 
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/ 9338899 3726897r^ 48350497,,, 4793, 

+ ni + C(3) h 

V 2099520 839808 1399680 ^^^ 58320 

/ 27670774337 14470577r^ 95617883401 

^ V 1414551600 ^ 765450 943034400 ' 

128,,^, 348701, \ 5397779543 265316737l7r4 



27^'' 17010 V 146966400 881798400 

35 9 6 877r\ ^ 17554601717,,^, 3655 ,,^, 84951877, 

In 2 + C(3) C(5) h 

229635 32659200 ^ ^ 10206^ ^ 2449440 

-'^h. (11) 

5103 ^ ' 

These results agree with the semi-analytical ones presented in Refs. |10|llj . 

As the authors of Refs. |10|llj . we conclude that the inclusion of the new four-loop 
corrections changes the actual values of the charm- and bottom-quark masses derived from 
QCD sum rules only very little, while it strongly decreases their theoretical uncertainties: 
approximately by factors of three and four in the case of charm and bottom, respectively. 



4 Conclusion 

We presented analytical results for the four-loop tadpole master integrals involving one 
non- vanishing mass that are called T54, Tq2, and T91 in the recent literature |10ll3j . for 
which only numerical values had been available so far. They represent essential ingredients 
for high-precision predictions of a number of observables of current phenomenological in- 
terest. As such an application, we considered the vacuum polarisation induced by a heavy 
quark in the four-loop approximation of QCD and expressed the first two coefficients of 
its Taylor expansion in the ratio of virtuality to heavy-quark mass in fully analytical 
form. The introduction of the characteristic combinations 64 and 65 of basic transcen- 
dental numbers turned out the be useful in order to compactify the expressions. Other 
applications include the four-loop QCD correction to the electroweak p parameter, which 
has just been presented in semi-analytic form [TU]. 
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